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IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

In re Patent Application of 

MYKLEBUST, H. et al. Atty. Ref.: 2810-17 

Serial No. uknown Group: 
US National Phase of PCT/NO00/00289 

Filed: March 7, 2002 Examiner: 

For: SYSTEM FOR PREDICTING THE OUTCOME OF AN IMAGINARY 
DEFIBRILLATOR SHOCK 



* * * 



March 7, 2002 

Assistant Commissioner for Patents 
Washington, DC 20231 

Sir: 

PRELIMINARY AMENDMENT 

In order to place the above-identified application in better condition for 
examination, please amend the application as follows: 
IN THE SPECIFICATION 

Please substitute the following paragraphs in the specification for corresponding 
paragraphs previously presented. A copy of the amended specification paragraphs 
showing current revisions is attached. 
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Page 1, before the first line, please insert as a separate paragraph: 

-This application is the US national phase of international application 

PCT/NOOO/00289 filed 6 September 2000, which designated the US.- 



IN THE CLAIMS 

Please substitute the following amended claims for corresponding claims 
previously presented. A copy of the amended claims showing current revisions is 
attached. 

4. (Amended) A system according to Claim 1 . w h e r e i n the algorithm for 
calculation the probability figure is a table look up in a m-dimensional table, where there 
for each table element is stored a numerical value for the probability figure, the table look 
up is determined from the value of a m-dimensional vector, the value of the tri- 
dimensional vector is diverted from the calculation of the energy of respective m 
different signal sequences that is represented on the output of m different digital filters, 
where the signal on the input of each digital filter is the segment of the ECG signal. 

6. (Amended) A system according to Claim 1 , w h e r e i n the calculation unit 
is connected to a data storage, the calculation unit is storing for each treatment 
parameters which describe the patient and parameters which describes the treatment, the 
calculation unit is connected to means for exchange of data, the exchange of data occur 
on a regular basis towards a central computer, where the calculation unit receives 
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optimised algorithm for calculation of the probability figure, and the computer receives 
information that is stored en the data storage. 

7. (Amended) A system according to Claim 1, w h e r e i n there is provided 
an optimised algorithm by first establishing of an updated set of empirical data consisting 
of information from a number of new patient treatments together with information from a 
number of earlier performed patient treatments, which all contain sequences of ECG 
where the outcome after shock are known; the optimised algorithm occur by iterative 
search after filter coefficients by m digital filters, where the filter coefficients are adjusted 
iterative in view of performance of a classification routine, where again the classification 
routine is adjusted iterative in view of performance and generality, where the 
performance is defined as the sum of sensitivity and specificity for classification of each 
of the ECG sequences to outcome classes ROSC and non-ROSC respectively, where the 
real outcome of shock is known for each ECG sequence, generality is fulfilled as the 
classification routine has the same performance for a arbitrary composite half of 
empirical material as for the rest of the empirical material, where measurement of 
generality and performance is provided in that each ECG sequence in the empirical 
material is expressed as a in-dimensional vector calculated from energy at the output of m 
digital filters, where the classification routine classifies each in-dimensional vector to one 
of the outcome classes ROSC, non-ROSC respectively, where the performance is 
measured as the sum of sensitivity and specificity of the classification routine, where a 
arbitrary composite half of empirical material has the same performance as the rest of the 
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empirical material, the optimised algorithm for calculation of the probability figure 
consist of a matrix having m matrix elements, where each matrix element express a 
probability figure, where the probability figure for each matrix element is provided by 
grouping ECG-sequences which is expressed by approximately identical m-dimensional 
vectors, where the occurrence of ECG which resulted in ROSC by shock plus the sum of 
occurrence of ECG which resulted in non-ROSC by shock constitute the probability 
figure for the matrix element, the in-dimensional matrix together with the filter 
coefficient constitute the optimised algorithm for calculation of the probability figure. 

8. (Amended) A system according to Claim 1, w h e r e i n the output of the 
analysis unit is connected to a receiver in the shape of a display unit. 

9. (Amended) A system according to Claim 1, w h e r e i n the receiver is a 
defibrillator. 

14. (Amended) A system according to Claim l,whereina device for 
indicating patient specific information and/or specific information regarding the 
treatment is connected to the analysis unit (2). 



-4- 



608845 



MYKLEBUST, H. et al. 
Serial No. uknown 

US National Phase of PCT7NO00/00289 



REMARKS 



1, 0 O 7 Q 5 s+S iii O 6 fj H- O E 



Attached hereto is a marked-up version of the changes made to the specification 
and claims by the current amendment. The attached page(s) is captioned " Version With 
Markings To Show Changes Made ." 

Respectfully submitted, 



MNL:ecb 

1 100 North Glebe Road, 8th Floor 
Arlington, VA 22201-4714 
Telephone: (703) 816-4000 
Facsimile: (703) 816-4100 



NIXON & VANDERHYE P C. 




Michelle N. Lester 



Reg. No. 32,331 
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VERSION WITH MARKINGS TO SHOW CHANGES MADE 

IN THE SPECIFICATION 

Page 1, before the first line, please insert as a separate paragraph: 
-This application is the US national phase of international application 
PCT/NO00/00289 filed 6 September 2000. which designated the US.-- 

IN THE CLAIMS 

4. (Amended) A system according to Claims 1-3 , w h e r e i n the algorithm 
for calculation the probability figure is a table look up in a m-dimensional table, where 
there for each table element is stored a numerical value for the probability figure, the 
table look up is determined from the value of a m-dimensional vector, the value of the Tri- 
dimensional vector is diverted from the calculation of the energy of respective m 
different signal sequences that is represented on the output of m different digital filters, 
where the signal on the input of each digital filter is the segment of the ECG signal 

6. (Amended) A system according to Claims l-£, w h e r e i n the calculation 
unit is connected to a data storage, the calculation unit is storing for each treatment 
parameters which describe the patient and parameters which describes the treatment, the 
calculation unit is connected to means for exchange of data, the exchange of data occur 
on a regular basis towards a central computer, where the calculation unit receives 
optimised algorithm for calculation of the probability figure, and the computer receives 
information that is stored en the data storage. 
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7. (Amended) A system according to Claims 1-6, wherein there is 
provided an optimised algorithm by first establishing of an updated set of empirical data 
consisting of information from a number of new patient treatments together with 
information from a number of earlier performed patient treatments, which all contain 
sequences of ECG where the outcome after shock are known; the optimised algorithm 
occur by iterative search after filter coefficients by m digital filters, where the filter 
coefficients are adjusted iterative in view of performance of a classification routine, 
where again the classification routine is adjusted iterative in view of performance and 
generality, where the performance is defined as the sum of sensitivity and specificity for 
classification of each of the ECG sequences to outcome classes ROSC and non-ROSC 
respectively, where the real outcome of shock is known for each ECG sequence, 
generality is fulfilled as the classification routine has the same performance for a 
arbitrary composite half of empirical material as for the rest of the empirical material, 
where measurement of generality and performance is provided in that each ECG 
sequence in the empirical material is expressed as a in-dimensional vector calculated 
from energy at the output of m digital filters, where the classification routine classifies 
each in-dimensional vector to one of the outcome classes ROSC, non-ROSC respectively, 
where the performance is measured as the sum of sensitivity and specificity of the 
classification routine, where a arbitrary composite half of empirical material has the same 
performance as the rest of the empirical material, the optimised algorithm for calculation 
of the probability figure consist of a matrix having m matrix elements, where each matrix 
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element express a probability figure, where the probability figure for each matrix element 
is provided by grouping ECG-sequences which is expressed by approximately identical 
m-dimensional vectors, where the occurrence of ECG which resulted in ROSC by shock 
plus the sum of occurrence of ECG which resulted in non-ROSC by shock constitute the 
probability figure for the matrix element, the in-dimensional matrix together with the 
filter coefficient constitute the optimised algorithm for calculation of the probability 
figure. 

8. (Amended) A system according to Claims 1-7, w h e r e i n the output of 
the analysis unit is connected to a receiver in the shape of a display unit. 

9. (Amended) A system according to Claims 1-8, w h e r e i n the receiver is a 
defibrillator. 

14. (Amended) A system according to Claims 1-+^, w h e r e i n a device for 
indicating patient specific information and/or specific information regarding the 
treatment is connected to the analysis unit (2). 
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SYSTEM FOR PREDICTING THE OUTCOME OF AN IMAGINARY DEFIBRILLATOR SHOCK 



The present invention regards a system for calculating the probability figure for the 
5 outcome of an immediate following defibrillator shock on the basis of properties of the 
heart measured during sudden cardiac arrest and resuscitation, as stated in the 
introduction to Claim 1 . 

Near 40% of all those who suffer sudden cardiac arrest will be able to survive if they 
10 receive good, lifesaving treatment immediately. When treatment is delayed, the chances 
of survival decrease, cf. the article by Holmberg S, Holmberg M: "National register of 
sudden cardiac arrest outside of hospitals " 1998 [1]. The treatment primarily consists of 
cardio-pulmonary resuscitation (CPR), which is administered until a defibrillator is in 
place. Thereafter, the treatment consists of alternating use of a defibrillator and CPR 
15 until resuscitation or until an ALS team arrives. (ALS = "Advanced Life Support") The 
latter also includes medication and securing of the respiratory passages as part of the 
treatment, cf. ILCOR, "Advisory statements of the International Liaison Committee on 
Resuscitation. Circulation" 1997;95:2172-2184 [6] 

20 Scientific papers in recent years point out a number of factors that affect the chances of 
survival: 

Time: The chance of surviving sudden cardiac arrest falls with time from 

heart failure until the first defibrillator shock is administered.[l] 
CPR: The chance of surviving increases when someone administers CPR 
before the defibrillator arrives.[l] 
Quality of Studies show that the quality of the CPR influences the survival- 
CP 111 (Cf. the publications by Wik L, Steen PA, Bircher NG. "Quality of 
bystander CPR influences outcome after prehospital cardiac arrest. 
Resuscitation" 1994;28:1 95-203[2]. 

Gallagher EJ, Lombardi G, Gennis P. "Effectiveness of bystander 
CPR and survival following out-of-hospital cardiac arrest". J Am 
Med Assoc 1995;274:1922-5[3] 

Van Hoyvegen RJ, Bossaert H. "Quality and efficiency of 
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bystander CPR. Resuscitation" 1993;26:47-52[4]) 

Timing of A study shows that when the duration of sudden cardiac arrest 
CPR and exceeds a number of minutes, the chance of survival will increase if 
defibrillator the ambulance personnel first administers a period of CPR before the 
treatment: defibrillator is used. (Cf Cobb L, et al. "Influence of cardiopulmonary 
resuscitation in patients with out-of hospital ventricular fibrillation. 
JAMA", April 7, 1999 - Vol 281, No 13 [5] 



In case of sudden cardiac arrest, the electrical activity in the heart (ECG) will indicate 
the state of the heart. Today's defibrillators measure and analyse ECG in order to 
classify the rhythm. If the rhythm is classified as Ventricle Tachyardia (VT) or 

5 Ventricle Fibrillation (VF), defibrillator treatment may have an effect. VT is often the 
precursor of VF. VF will as time goes by and the energy and oxygen reserves of the 
heart muscle are depleted tend towards Asystole, a rhythm characterised by very little or 
no electrical activity. The purpose of the defibrillator treatment is to restore the 
organised electrical activity of the heart and the associated blood pressure and blood 

10 circulation. This is often denoted ROSC - "Return of Spontaneous Circulation", and is 
the first step towards survival. 

Only a fraction of the shocks delivered actually result in ROSC. Most shocks today do 
not give ROSC , cf. the publications Gliner BE et al. "Treatment of out-of hospital 

15 cardiac arrest with a Low-Energy Impedance-Compensating Biphasic Waveform 
Automated External Defibrillator" [7], Sunde K, Eftestol T, Askenberg C, Steen PA. 
"Quality evaluation of defibrillation and ALS using the registration module from the 
defibrillator. Resuscitation" 1999 [14]. In general, it can be said that the chance of 
ROSC is at its greatest immediately after sudden cardiac arrest, when the heart muscle 

20 still possesses energy reserves and oxygen. Many patients achieve ROSC after 

alternating use of shocks and CPR. The disadvantages of having to give many shocks 
are several: First of all, no CPR will be given during the shock treatment, a factor that 
further aggravates the situation for the vital organs, particularly for the brain. 
Furthermore, it has been shown that the heart muscle is also damaged by the shocks. 
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and that the damage increases with the number of shocks and the amount of energy, cf. 
the publications Ewy GA, Taren D Bangert J et al. "Comparison of myocardial damage 
from defibrillator discharges at various dosages. Medical instrumentation" 1 980; 14:9- 
12. [16]. For the patient, the ideal would be to be given only one shock, and for this 
5 shock to give ROSC. 

Thus, for many patients, it is crucial that the administration of CPR be effective, so as to 
revitalise the heart through supplying a flow of blood through the heart muscle, cf. the 
publication Michael JR et al. "Mechanism by which augments cerebral and myocardial 

10 perfusion during cardiopulmonary resuscitation in dogs. Circulation" 1984;69:822-835. 
[17]. This revitalisation can be indicated through ECG measurements, where ECG 
characteristics such as form, spectral flatness measurements, frequency, amplitude, 
energy etc. is seen to change back towards their initial values such as they would have 
been immediately after the heart action was suspended, cf. the publications Eftest0l T, 

15 Aase, SO, Husoy JH. "Spectral flatness measure for characterising changes in cardiac 
arrhythmias." Computers in Cardiology, [15] and Noc M, Weil MH, Gazmuri SS, 
Biscera J and Tang W. "Ventricular fibrillation voltage as a monitor of the effectiveness 
of cardiopulmonary resuscitation" in J Lab Clin Med, September 1994 [13]. This 
revitalisation will increase the probability of the next shock giving ROSC. 

20 

Unfortunately, not everyone survives. For many, the reason behind the sudden cardiac 
arrest is such that resuscitation is impossible. Furthermore, the time factor and the 
quality of the treatment will also play a part and affect the survival. 

25 Today's resuscitation guidelines describe a protocol that is the same for everyone, 
regardless of sex. race, how long the heart action has been suspended, whether a 
member of the public has given CPR etc. The means of resuscitation are primarily CPR 
and defibrillator treatment, and later also medication administered by lifesavers who 
have been given special training in this area. Today's protocol is such that if the first 

30 three shocks have no effect, CPR is to be given for 1 minute, then three more shocks, 
and so on. As it takes about one minute to give three shocks, the patient will be without 
CPR for half of the time. 
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Literature and other patent applications describe technology, the object of which is to 
guide the lifesaver in the choice between CPR and defibrillator treatment. Brown et al in 
US patents no. 5 683 424 and 5 571 142 [10] describe a system that, based on spectral 

5 measures in VF, instructs the lifesaver to either give CPR or give a shock. A separate 
analysis of this method, where the method has been tested on human VF, yields results 
that show the method to have a low specificity, i.e. that the method will only to a 
limited degree reduce the number of unnecessary shocks. Noc M, Weil MH, Tang W, 
Sun S, Pernat A, Bisera J. "Electrocardiographic prediction of the success of cardiac 

10 resuscitation" in Crit Care Med 1999 Vol 27 No 4 [12] describe a similar system, based 
on an animal model, which links the mean amplitude and dominant frequency of VF to 
the outcome of the defibrillator shock. Both of these methods aim to advise against 
defibrillator use as long as the condition of the heart is such that a shock is assumed not 
to have an effect, and instead use CPR. Both methods define absolute criteria based on a 

is limited number of observations from a defined group of patients or animals. 

The object of the present invention is to seek to constantly optimise the treatment 
through: 

20 a) By means of equipment connected to the patient, measuring the properties of the 

heart during sudden cardiac arrest, and, based on relevant treatment, knowledge of 
the patient and knowledge of comparable conditions and outcomes of treatment 
carried out previously, calculate the probability figure of achieving ROSC - Prose. 

b) Presenting the probability figure or using the probability figure in support of the 
25 decision on further treatment. 

c) Relaying registrations for each treatment, along with the outcome of each shock, 
to a centrally located computer, and using this experience to improve the basis for 
Prose calculations, so as to improve the confidence of the next calculation. 

d) Possibly registering the measured CPR parameters and, together with 
30 information regarding medication given, looking at the development of the 

probability figure for the purpose of identifying effective lifesaving. 
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e) Possibly providing feedback to the user regarding which CPR parameters and 
what medications have been identified as effective, or alternatively through using 
this information in order to instruct the user in effective CPR. 

5 To the extent the word "possibly" is used, this is done in order to take into consideration 
the fact that there are many categories of users and lifesavers, and based on their 
professional level and economy, it is appropriate to use the solutions that best serve the 
purpose. 

10 The object of the invention is to contribute towards giving the patient a treatment that is 
better suited to the individual, and which gives a greater chance of survival. The use of 
empirical data makes it possible to allow for differences and the ever-changing patient 
groups and treatment. To the extent that not all shocks give ROSC, the system attempts 
to take into consideration factors that affect or may affect the probability of ROSC: 
15 • Properties of the heart that may be observed via the electrodes, and which 

express the metabolic condition and the pumping action. These properties 
develop with time during sudden cardiac arrest, but may be partially reversed 
through the use of medication and CPR. 

• Type of defibrillator shock and energy selection produce varying effectiveness. 

20 • Physical conditions. A big patient will receive a lower current density through the 
heart than a small patient, for the same energy selection. 

• Patient information. The system takes into account the fact that there may be a 
difference between men and women, based on the fact that over 70% of those 
who suffer sudden cardiac arrest are men. In certain parts of the world, the 

25 duration of life has increased, so that the number of elderly who suffer sudden 

cardiac arrest is increasing. These may very well have received treatment over a 
period of time, both medications, surgical procedures and aids such as 
pacemakers, all of which may have an effect on the Prose. 

• Geography, race. The system further takes into consideration the fact that there 
30 may be differences based on lifestyle and genetic conditions, just as life 

expectancy varies greatly depending on geography and race. 
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The above is provided by means of a system of the type mentioned initially, the 
characteristics of which are stated in Claim 1. Further characteristics of the invention 
appear from the remaining, dependent claims 

Using Prose to optimise the treatment may be done in several ways. For advanced users, 
the most expedient will be to present the parameter graphically versus time, as a trend 
curve. This will immediately provide a direct indication of the state of the heart, and 
also indicate the effect of medication and CPR. For groups who are not trained in 
relating to this type of information, the most appropriate thing will be to provide 
automatic decision support in the question of whether or not to give CPR, in which way 
CPR should be given, or whether shocks should be given. The principle of simple 
decision support may be: 

• If Prose is less than a limit, Grose, CPR is recommended. Otherwise, a number 
of sequential defibrillator shocks are recommended. 

• CPR is recommended until positive changes in the Prose level out, but no further 
than T minutes before a number of defibrillator shocks are recommended. 

The following will describe the invention in greater detail, with reference to the 
drawings, in which: 

Fig. 1 shows system components consisting of one, alternatively several, 

computers in a network that communicates with a number of positioned 
analysis units. 

Fig. 2 shows the block diagram for a defibrillator with a built-in analysis unit. 

Fig. 3 shows the elementary flow diagram for information. 

Fig. 4 shows an apparatus with electrodes connected to the patient's chest, in 

positions on the chest that are normally used for delivery of a 

defibrillator shock, as well as for measuring ECG in accordance with 

standard derivation II. 
Fig. 5 shows a flow diagram for development of coefficients to optimal filter. 
Fig.6 shows a flow diagram for development of a classification that fulfil the 

requirement of generality. 
Fig. 7 shows a general block diagram of the invention with focus on the 

analysis unit. 
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The system consists of one, alternatively several, computer(s) 1 in a network that can 
communicate with a number of positioned analysis units 2. These may either be 
integrated into equipment (Ul, U2..) such as defibrillators or ECG monitors, or they 
may occur in or as a support product used during the resuscitation attempt The analysis 
5 units 2 generally operate independently of the computers 1 , however after use, the 
analysis units will deliver field data to the computer 1, and possibly receive adjusted 
algorithms for calculation of property vector and/or Prose 

The analysis unit 2 normally has subsystems associated with it, cf. Fig. 2: 
10 Some of these subsystems are standard in equipment such as defibrillators and ECG 
monitors, and these are as follows: 

Reference number 3 denotes a system for measuring ECG, connected to electrodes E. 
For measuring and digitalize of bio-electrical signals is used said electrodes connected 
to skin of the patient in accordance with standard derivation II, cf. Fig. 4, which shows 

15 the device with electrodes connected to the patient's chest in locations that are normally 
used for delivering defibrillator shocks, as well as for measuring ECG in accordance 
with standard derivation II. In this derivation, ECG will essentially express the electrical 
activity in the longitudinal axis of the heart. Reference number 4 denotes a system for 
checking the electrode connection to the patient, reference number 5 denotes a device 

20 for high voltage and shock delivery, if integrated into a defibrillator, reference number 7 
denotes an algorithm for classification of ECG, if integrated into a defibrillator, 
reference number 7 denotes a microprocessor system and software, reference number 8 
denotes computer memory, reference number 9 denotes a user interface, reference 
number 10 denotes an energy supply, reference number 1 1 denotes a system for delivery 

25 and receipt of data from external equipment. 

Subsystems 3-1 1 are standard equipment in defibrillators and ECG monitors, and will 
therefore not be described further in this specification. 

30 Analysis unit 2 consists of the following units: 
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Unit 12 for determining one or more properties of the heart that is combined in a 
property vector, and based on this calculate the probability of ROSC, Prose, for the 
patient who is connected up. Possibly also a module 13 for determining the blood flow 
through the heart, based on the measured impedance and the change of the impedance 

5 between the electrodes as a function of the pumping action of the heart and the 
expansion of the lungs, a possible module 14 for registering CPR parameters from 
sensors S, a possible module 15 for registering patient specific information, a possible 
module 16 for registering any medication administered; and a possible module 17 for 
correlating positive changes in Prose with information regarding the treatment given, 

10 and display or use this information to guide the treatment. 

The analysis unit may be realised as a self-contained microprocessor unit, or it may be 
realised through the microprocessor unit in the equipment. 

15 Each computer consists of the following subsystems: 

(a) Hardware, (b) operating system, (c) software and interface for communication in a 
network (d) database for field data, (e) algorithm for calculation of property vector, 
including a classification system for predict the outcome of an immediate following 
20 defibrillator shock, (f) system for calculation of Prose, (g) algorithm for correlating 
changes in Prose with information regarding patient and treatment, and (h) system for 
delivery and receipt of data from positioned defibrillators. 

With regard to the computers, the subsystems of hardware, operating systems, software 
25 and interface are of a generic nature, and will not be described in greater detail. 

A more detailed description of the system, first computer 1. 

(d) The database for field data consists of a large amount of patient episode data, and 
30 contains: 

- Patient information: Sex, age, weight, race etc. 

- Geographical information 
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- Information regarding each defibrillator shock: Curve shape, energy, timing 
versus VF. 

- For each shock: 

- Preshock ECG 

5 - Preshock CPR data 

- Preshock medication data 

- Preshock impedance data 

- Postshock ECG 

- Postshock impedance data 

10 - Annotation of ROSC/Non-ROSC, with outcome rhythm for each shock 

(e) The algorithm for calculation of the property vector (v) makes use of mathematical 
methods in order to characterise the condition of the heart based on a recording of a bio- 
medical signal (x). The bio-medical signal is preferably ECG, but also other type of 
is signals such as signals derived from the change of the impedance, sound or pattern of 
the movement may be considered. 

The algorithm for calculation of the property vector is hereafter denoted as v(x). v(x) 
which is used on ECG from empirical data provide us two set of property vectors: 
20 An set, VI, containing nl property vectors where the outcome of the shock is ROSC, 
and an set, V2, containing n2 property vectors where the outcome of the shock is none- 
ROSC. 

In general, v(x) is defined as an operator that operates on an ECG sequence, x, 
25 consisting of TV samples, which generates a property vector, v, consisting of M vector 
elements that ideally takes care of the information in jc that separates the group of x that 
results in ROSC, XI, from the group of jc that results in non-ROSC, X2, The methods of 
property extraction are innumerable, and the literature describes some of these, which 
can be roughly divided into time and transform domain methods, where the object is to 
30 structure x in a manner that is appropriate for property extraction. Among preferred 
methods are: 
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Optimised digital filters determined by L filter parameters that divides x into M 
channels. The energy from each of these channels is calculated, so as to make 
the property vector consist of M elements. These types of filters are described 
inter alia by T. Randen "Filter and Filter Bank Design for Image Texture 
Recognition" in a thesis of NTNU, Oktober 1997 where the filters are optimised 
in order to achieve the best possible recognition of the different textures. For 
present purpose is the optimised filters found by using a numerical gradient 
search algorithm (T. Coleman, M. A. Branch and A. Grace, , Optimization 
Toolbox for Use with MATLAB, The Math Works Inc, 1999) to achieve the best 
possible separation of the ROSC group from the non-ROSC group. Separation 
ability is measured by the sum of sensitivity (degree of correct recognition of 
ROSC) and specified (degree of correct recognition of nom-ROSC) By a given 
iteration in optimisation is this performance measured, and the set of parameters, 
which define the filters, is adjusted in the direction corresponding the increase in 
performance. This procedure is repeated until maximum performance is reached. 
Spectral measures that are calculated on the basis of the estimate of the power 
density spectrum (PSD) of x. The PSD can be estimated through use of Fourier 
transforms. Based on the PSD, characteristics are calculated that express the 
frequencies at the centre and the maximum points of the PSD. In addition, the 
flatness and energy of the PSD is also characterised 

Examples of other methods for property extraction are: 

- wavelet analysis, 

- neural networks 

The relation between VI and XI, V2 and X2 respectively are as following: XI 
containing an set of nl ECG sequences, which, when used on v(x), provides a number 
of property vectors VI, which all is belonging to the outcome class ROSC (wl). X2 
containing an set of n2 ECG sequences, which, when used on v(x), provides a number 
of property vectors VI, which all is belonging to the outcome class non-ROSC (w2). 
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(f) A system for calculation of the Prose function is based on pattern recognition theory, 
and forms the second element of the classification system. In this context, the term 
classes is defined as the collection of measurements of the condition of the heart that 
corresponds to 

- ROSC (vW) 

- non-ROSC (w2) 

The property vectors of the two classes are statistically described by 

- P(wi), i=l 9 2 , which is the a priori probability of the two classes. I.e. 
before a measurement is made, the probability of one or the other 
outcome is known through the respective a priori probabilities. 

- p(v\wi) are the class specific probability density functions. These express 
how the measurements within the given classes are distributed. p(v) 
expresses the compound probability density function for the 
measurements, and is given by adding up the class specific probability 
density functions weighted by the associated a posteriori probabilities. 

- P(wi\v) are the a posteriori probability functions for the two classes. 
These functions express the probability of a given measurement 
belonging to wi. Bayes formula expresses P(wi\v) as a function of the 
above probability functions . 

- P(wi|v)=P(wi)*p(v|wi)/(p(wl)+P(w2)*p(v|w2)) 

- The sum of the a posteriori probabilities for a given v is always 1 . 

In the case of a given measurement, v, one wishes to determine the class allocation wl 
or w2. It has been proven that the expected probability of misclassification is minimised 
by selecting the wi that corresponds to the maximum P(wi\v). It is further possible to 
define (make an estimated choice of) the cost of all types of misclassification, such that 
the expected risk of a given misclassification is given by the product of the cost and the 
a posteriori probability of the true class. The expected risk of misclassification can then 
be minimised by classification is a class corresponding the product with smallest value 
is selected. 
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In most cases, the statistics of the property vector are not known. These quantities must 
then be estimated before Prosc(v) can be produced. The pattern recognition theory 
describes a multitude of methods for this, which are based on measurements (practice 
data) that are examples from the various wi. Some examples: 
5 - Histogram techniques, which divides the outcome space into hypercubes 

in which the probabilities within each of these are calculated on the basis 

of the number of occurrences of the different classes within the given 

hypercube. This corresponds to the method used herein. In the following 

is described how statistic quantity is estimated. 
10 We will start defining the quantities: 

n^total number of observations in the empirical material. 

nl= total number of observations corresponding ROSC outcome. 

n2 = total number of observations corresponding non-ROSC outcome. 

njl== total number of observations corresponding ROSC outcome within 
!5 hypercube no. j. 

n j2= total number of observations corresponding non-ROSC outcome 

within hypercube no. j 

We have n=nl+n2. Estimate for a priori probability will then be 
20 A P(wi)=ni/n, i=l,2. 

The local estimates (within hypercube j) for the class specific probability 
function will then be 
A p(v|wi)=nji/ni, 1=1,2. 

The local estimates for a posteriori probabilities is calculated in respect 
2 5 of the Bayes formula inserted estimate for a priori probability and the 

local class specific probability density functions. Se R. J. Schalkoff. 

Pattern recognition: Statistical, structural and neural approaches. John 

Wiley & sons, New York (NY), 1992 

A P(wi|v)=nji/(njl+nj2), i=l,2 
30 - Radial base functions, in which the probabilities at a given point are 

calculated on the basis of the contribution from surrounding practice data 

from the different classes. The contributions decrease with distance. 
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- Parametric modelling, in which a mean value and dispersion for the 
different classes are used to produce analytical probability models. 

- Neural networks, learning vector quantization and nearest neighbour 
classification are some other central methods within the pattern 

5 recognition theory. 

It is important that a given classifier be tested on a set of observations (test set) 
independently of the practice data (practice set), in order to check that the classifier 
yields the expected results. The demand is that there is consistency between practice and 

10 testing, that the classifier fulfils the requirement of generality. With generality means: 
By dividing the empirical data in two parts and letting the one part represent a set of 
data called practice set and the other part represent a set of data called test set, the 
generality is defined as following: Decision limits used on all of the property vectors in 
each set of data for classification of the outcome, which provides approximately same 

15 performance (the sum of sensitivity and specificity) for both set of data fulfil the 
requirement of generality. These decision limits occur through an iterative process 
where the practice set is included in the calculation of the decision limit, see fig. 6. 

Those measurements v that correspond to the ROSC outcome belong in wl. The 
20 probability of a given measurement, v, belonging in wl is given by P(wl\v). In other 
words, this probability function expresses the probability figure Prose of ROSC for a 
given measurement v. 

Prosc(v)=P(wl\v) 

25 

As mentioned previously, different property vectors, v, can be calculated by means of a 
countless number of methods. Which methods and which dimension, M, is suitable for 
expressing Prosc(v) is assessed on the basis of the expected risk in the case of 
misclassification for each method. The method that minimises this risk is the most 
30 appropriate for expressing Prosc(v). Fig.6 shows a flow diagram for an iterative 

development of algorithm for calculation of the property vector v. Basis for the iterative 
development is empirical data. As the amount of empirical data increase is this iterative 
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process repeated so that the ability of the property vector to predict outcome classes is 
increased. The iterative adjustment of the decision limit is also included so that the 
requirement of generality is fulfilled. 

5 (g) The algorithm for correlating changes in Prose with information regarding the 

patient and the treatment is mainly for scientific purposes. The defibrillator to guide the 
user during lifesaving may later use the results from the correlation. 

Prosc(v) has been provided as described under points (d) and (e). In this analysis, ECG 
10 sequences are extracted from the patient material, so that the ECG sequences describe a 
course of treatment that is as uniform as possible, with the complete patient material 
seen under one. Examples of such a course of treatment may be 

- CPR sequences 

- "Hands off intervals, for instance a period for defibrillator rhythm 
15 analysis up to the shock, after a CPR period. 

In these ECG sequences, corresponding Prosc(v) sequences are calculated as described 
under points (d) and (e). Consequently, the change in Prosc(v) , DProsc, is calculated 
for each sequence. DProsc is grouped on the basis of those treatment characteristics that 
20 are of interest with regard to the effect of the treatment. As an example, one can group 
DProsc with regard to the following treatment characteristics, singly or in combination: 

- Different compression frequencies, compression depths, duration of chest 
compression 

Degree of ventilation 
25 - Medication 

- Physiological measurements such as blood flow measurements, blood 
pressure etc. 

Where significant differences in DProsc occur for dissimilar treatment conditions, this 
information may be used to identify advantageous treatment methods. This information 
30 may be utilised through the person giving the treatment being given feedback regarding 
good and poor treatment. 
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(h) A system for delivering and receiving data from positioned analysis units. Here, no 
special requirements apply. The exchange of data can take place directly through use of 
memory modules such as PCMCIA, cordlessly by means of IR or RF communication, 
via networks such as the Internet, or by a direct connection between communication 
5 ports in the equipment and the computer. The most practical method these days is to 
have the analysis unit 2 communicate directly with the computer 1 via a local PC that it 
can communicate with, and to have the local computer pass the data on via the Internet. 

Detailed description of the analysis unit 2. 

10 

System 12, algorithm v(x) for calculating the property vector (v) and algorithm for 
calculating the probability of ROSC, Prose, for ECG from the patient who is connected 
up. 

- v(x)is a set with calculation, which combined form a property vector. The 
15 calculation is a set of energy calculation within a determined frequency 

band (optimised filter) or a set of parameters diverted from effect density 
spectrum or a combination of this. 

- The algorithm for calculation of Prose is typical a matrix where the 
number of dimension is consistent with the number of dimensions within 

20 the property vector, where each matrix element is containing a numerical 

value for pROSC. The numerical value is downloaded from the computer 
1 . The matrix is provided in that there for each ECG sequens within the 
empirical material is calculated a property vector. Property vectors that 
are approximately identical are grouped thereafter together in a number 

25 of sets, where each set is assigned a matrix element. For an set, j, the 

occurrences of ROSC is counted, njl, towards the occurrences of non- 
ROSC, nj2. The ratio between the numbers of {ROSC} towards the 
number of {ROSC+non-ROSC} define the numerical value of pROSC 
for the set j which is njl/(njH-nj2). This correspond the local estimates 

30 for a posteriori probability for the class wl as described by use of 

histogram technique. The numerical value for each set is then stored in 
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that matrix element that correspond to the elements of the property 
vector. 

The module 13 for calculating blood flow through the heart based on the measured 
5 impedance and the change of the impedance between the electrodes as a function of the 
pumping action of the heart and the expansion of the lungs: 

- The value of the measured impedance, Zo, measured by means of an 
approximately constant alternating current, informs the analysis unit 1 of 
the impedance between the electrodes, and can be used to replace system 

10 4. 

The impedance change between the electrodes will be proportional to the 
change in the set of air in the lungs plus the working volume of the heart. 
The change due to air dominates. By looking at the signal between two 
ventilations, or by first filtering out the ventilation, it will be possible to 

15 estimate the working volume on the basis of the formula 



AV= AZpL 2 
Zo 2 

20 

This formula is universally known, and is used within Impedance Cardiography. 
A7. is the impedance change, p is the resistivity of the blood, L is the distance between 
the electrodes, and Zo is the numerical value of the impedance. A simplification of this 
formula is preferable: 

25 AV= 

Zo 2 

Here, A: is a constant. This measurement will indicate to what degree the blood is 
flowing, and will contribute towards characterising the condition of the heart in VF/VT, 
30 and will furthermore indicate ROSC in case of a successful defibrillator shock. 

Module 14 for measuring and registering CPR parameters. 
Relevant CPR parameters are: 
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- Inflation time and inflation volume are measured by looking at the 
impedance change between the electrodes. This change is several times 
greater than the change that takes place as a function of the blood stream 
from the heart, and is proportional to the amount of air in the lungs. The 
principle is known from other diagnostic equipment. 

- Compression rate and compression proportion (the ratio between 
compression time and relaxation time) during chest compression may be 
measured by looking at the impedance change between the electrodes, or 
by placing a sensor at the compression point on the patient's chest. This 
sensor can contain a pressure switch, or a dynamometer or an 
accelerometer. 

- Compression depth is calculated on the basis of signals from an 
accelerometer placed at the compression point. 

- Time between inflation and chest compression, and time between chest 
compression and inflation. 

- Proportion of CPR relative to the total treatment time 

- Amount of compression, the sum of the product between the duration and 
depth of the compression. 

Module 15 for indicating patient specific information. 

This information can be passed to the analysis unit 1 e.g. by dedicated push buttons, or 
it may come in from an external source such as a patient database or a patient journal on 
a PC / handheld computer. Relevant information is: 
Geographical area 

- Age 

- Sex 

- Weight 

- Race 



Module 16 for indicating medication and dosage given. 
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This information can be passed on to the analysis unit 1 through dedicated push buttons, 
from a patient journal on a PC or from other devices that log the use of medication. 
Relevant medicines are 

- Epinephrine 
5 - Lidocaine 

Bretylium 

- Magnesium sulphate 

- Procainamide 

- Vasopressins 

io - Thrombolysis medication 

- and so on. 

Module 17 for correlating positive changes in Prose with information regarding the 
treatment given, and displaying or using this information to guide the treatment. 
15 - The system identifies and registers a period of Prose with positive 

change. At the same time, the system identifies and registers the average 
of each CPR parameter measured for a period of time prior to the change 
and during the change, and if applicable, what medication was given 
during the same period. 
20 - If a new period of improved positive change occurs, the identification of 

CPR parameters and medication is repeated. 

- This information can be displayed on the defibrillator screen, or it may be 
used to produce voice messages that guide the user to deliver CPR with 
parameters that are in accordance with the registration. 

25 - This information will also be of great importance to research, with a view 

to optimising the guidelines for CPR treatment and training. 

To summarise, the invention can be described with reference to Figure 3, which 
schematically shows the flow of information between the central computer 1 and the 
30 analysis unit 2 in the positioned equipment in an application of the invention. 



The computer 1 contains empirical data from previous resuscitation attempts, where the 
outcome of the resuscitation attempt is known. The main ingredient in the empirical 
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basis is the ECG and the associated outcome after a shock (ROSC/non-ROSC). 
Additional empirical data impart nuances to the relationship between outcome, 
treatment and patient specific factors. 

5 An algorithm v(x) for calculating of a property vector v is developed based on empirical 
material. This algorithm is adapted so that the ability of the property vector to predict 
the outcome (ROSC/non-ROSC) is optimised within the demand of generality. 

A property vector v is thus characterised by that it represents ECG segments that is 
10 succeeded by ROSC after shock characteristically different from ECG segments that is 
succeeded by non-ROSC after shock. Such property vector is provided through an 
iterative process where the parameters for each iterative step is being adjusted until the 
best performances are reached, see fig. 5. 

is By grouping the property vectors in the empirical data set that are approximately 

identical and adding up the number of cases of ROSC and the number of cases of non- 
ROSC for each of these sets, then the probability figure Prose for ROSC is expressed 
statistically relative to the value of an in-coming property vector which represents the 
ECG segment Prose is desired to be calculated for. 

20 

The property vector of a ECG segment consist of a number of element which together 
characterize the ECG segment with a view to the outcome after a shock. Let the 
property vector consist of (n+p) elements that together describes one determined ECG 
segment. Then the property vector will be n energy calculations of the output signal of n 

25 different digital filters, where the ECG segment constitute the input signal plus p 
parameters that express characterizing features by the power density spectrum of the 
ECG segment. Preferred characterizing features is: Frequency by the centre of gravity, 
frequency by the maximum point, spectral flatness measurements, and spectral energy. 
These parameters will to a limited degree overlap each other in view of information 

30 contents. In order to reduce the set of data without to reduce the number of information 
may for example be used principal component analysis. By such a method may a big set 
of data, which is described of (n+p) elements in each vector be reduced to one set of 
data where each vector only have m elements; m<(n+p). 

35 A practical way of expressing this statistical interrelationship is through a Prose 
function, which is a substitute for all the empirical data, but which mathematically 
expresses the same relationship between the property vector and Prose. 
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This function is entered into the programme code of the analysis unit, so that when this 
receives a segment of ECG, the analysis unit will first perform the same calculation of 
the property vector as that performed by the computer, and then use the property vector 
5 as input to the Prose function in order to calculate the probability figure of an immediate 
following defibrillator shock giving ROSC. 

The analysis unit may furthermore impart nuances to, and thereby provide a more 
accurate probability figure, by also utilising knowledge of the treatment and the patient, 
10 seen in relation to empirical data. 



The ever-changing forms of treatment and patient characteristics necessitate a 
continuous update of the empirical basis. This is achieved by each analysis unit passing 
on its experience regarding each shock to the central computer, where the central 
is computer repeats the grouping of the property vectors, recalculates the Prose function 
and passes the result back to the analysis unit. 
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Claims 
1. 

A system for calculating probability figures for the outcome of an immediate following 
5 defibrillator shock resulting in return of spontaneous circulation (ROSC) 
wherein 

an analysis unit is connected with a module, which is measuring bio-electrical signals 
from electrodes connected to a patient, 

the analysis unit is provided to organise the bio-electrical signals continuous into 
10 segments, 

the analysis unit is provided for each segments to calculate a combination parameter 
that characterise the condition of the heart 

the analysis unit is provided to by means of comparing for each combination of 
parameters to find a corresponding combination of parameters from earlier made 
15 defibrillator treatments, where there for each combination of parameters is assigned a 
probability figure, where the probability figure expresses the number of defibrillator 
shocks that results in ROSC relative the total number of defibrillator shock for each 
combination of parameters, and 

the analysis unit has an output for the probability figure, 

20 

2. 

A system according to Claim 1 , wherein 
the bio-electrical signals is ECG signals. 

25 3. 

A system according to Claim . |-1 w h e r e i n 

the analysis unit (2) is provided to calculate the probability figure by means of an 
algorithm. 

30 4. 

A system according to Claims 1-3 , wherein the algorithm for 
calculation the probability figure is a table look up in a m-dimensional table, where 
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there for each table element is stored a numerical value for the probability figure, the 
table look up is determined from the value of a m-dimensional vector, the value of the 
m-dimensional vector is diverted from the calculation of the energy of respective m 
different signal sequences that is represented on the output of m different digital filters, 
where the signal on the input of each digital filter is the segment of the ECG signal. 

5. 

A system according to Claims 4, wherein 

the value of the m-dimensional vector is diverted from calculation of flatness, energy, 
frequency by the centre of gravity and frequency by the maximum point of a power 
density spectrum, where the power density spectrum is diverted from the ECG signal 
segment. 

6. 

A system according to Claims 1-5, wherein the calculation unit is 
connected to a data storage, the calculation unit is storing for each treatment parameters 
which describe the patient and parameters which describes the treatment, the calculation 
unit is connected to means for exchange of data, the exchange of data occur on a regular 
basis towards a central computer, where the calculation unit receives optimised 
algorithm for calculation of the probability figure, and the computer receives 
information that is stored en the data storage. 

7. 

A system according to Claims 1-6, wherein 

there is provided an optimised algorithm by first establishing of an updated set of 
empirical data consisting of information from a number of new patient treatments 
together with information from a number of earlier performed patient treatments, which 
all contain sequences of ECG where the outcome after shock are known; the optimised 
algorithm occur by iterative search after filter coefficients by m digital filters, where the 
filter coefficients are adjusted iterative in view of performance of a classification 
routine, where again the classification routine is adjusted iterative in view of 
performance and generality, where the performance is defined as the sum of sensitivity 
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and specificity for classification of each of the ECG sequences to outcome classes 
ROSC and non-ROSC respectively, where the real outcome of shock is known for each 
ECG sequence, generality is fulfilled as the classification routine has the same 
performance for a arbitrary composite half of empirical material as for the rest of the 
empirical material, where measurement of generality and performance is provided in 
that each ECG sequence in the empirical material is expressed as a m-dimensional 
vector calculated from energy at the output of m digital filters, where the classification 
routine classifies each m-dimensional vector to one of the outcome classes ROSC, non- 
ROSC respectively, where the performance is measured as the sum of sensitivity and 
specificity of the classification routine, where a arbitrary composite half of empirical 
material has the same performance as the rest of the empirical material, the optimised 
algorithm for calculation of the probability figure consist of a matrix having m matrix 
elements, where each matrix element express a probability figure, where the probability 
figure for each matrix element is provided by grouping ECG-sequences which is 
expressed by approximately identical m-dimensional vectors, where the occurrence of 
ECG which resulted in ROSC by shock plus the sum of occurrence of ECG which 
resulted in non-ROSC by shock constitute the probability figure for the matrix element, 
the m-dimensional matrix together with the filter coefficient constitute the optimised 
algorithm for calculation of the probability figure. 

8 

A system according to Claims 1 wherein 

the output of the analysis unit is connected to a receiver in the shape of a display unit. 
9. 

A system according to Claim 8, wherein 
the receiver is a defibrillator. 

10. 

A system according to Claims 1, wherein the receiver of the probability 
figure is an algorithm for decision support for the choice of treatment. 
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11. 



A system according to Claims 1, wherein 

the analysis unit (2) identifies periods of positive change in the probability figure 
together with parameters that characterise the treatment, and passes on the numerical 
5 value of the positive change in the probability figure, together with the mean value of 
each treatment parameter over the period, to a receiver. 



A system according to Claims 1 , wherein 
10 the receiver of the numerical value of the positive change in the probability figure, 
together with the mean value of each treatment parameter over the period, is a display 
unit. 



15 A system according to Claims 1, wherein 

the receiver of the numerical value of the positive change in the probability figure, 
together with the mean value of each treatment parameter over the period, is an 
algorithm for decision support for the choice of treatment. 

20 14. 

A system according to Claims 1-13, wherein 

a device for indicating patient specific information and/or specific information 

regarding the treatment is connected to the analysis unit (2). 
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